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Self-organization effects related to light amplification in the collective atomic recoil laser system 
with the driven atoms confined in a harmonic trap are investigated further. In the dispersive 
parametric region, our study reveals that the spontaneously formed structures in the phase space 
contributes an important role to the light amplification of the probe field under the atomic motion 
being modified by the trap. 
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I. INTRODUCTION AND SIMULATION RESULTS 

Optical amplification in a strongly driven sodium atoms system without population inversions has been observed 
in the experiment [l| where the atomic recoil effect and Doppler shift were purposely avoided by a specific systematic 
configuration. When the pumping field is off-detuned from the atomic transition frequency loq, the experimental 
pump-probe gain spectrum exhibits a normal absorption-gain profile similar to that of Fig. 1(a) which is drawn from 
a direct simulation and has been obtained by Bonifacio et al in their original work The vertical ordinate of 
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FIG. 1: Gain profile G(t, A21) of the probe field versus pump-probe detuning A21 with parameters v = 0, Y = 1, k = 0.01, 
p — 3, A2 = 2, A20 = — 15 and r = 5. (a) the atomic motion is totally removed from the model; (b) the motion effect is 
included (traditional CARL). 



Fig.l stands for a relative gain of the probe field and the horizontal ordinate A21 denotes the detuning between the 
pump and the probe field. However, this anti-symmetry gain profile around A21 = can be drastically changed 
(Fig. 1(b)) by the atomic recoil motion which plays an important role in the amplification mechanism of the collective 
atomic recoil laser (CARL) 0. As shown in Fig. 1(b), the gain peak within the traditional Madey region 3] is greatly 
enhanced, and which has been a signature of light amplification induced by cooperative atomic recoil effect. Being 
modified by the recoil motion and exchanging energies with one common pump field, the atoms will exhibit some 
kinds of self-organization behavior Q or dynamic synchronization [f| which induces phase transitions in the relevant 
system @|3- Indeed, in the traditional CARL system, the atomic density grating as well as the population inversion 
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grating on the scale of the electromagnetic wavelength emerges spontaneously during the light amplification process 
ra. However, some theoretical simplicities of the traditional CARL model, as have been pointed by some authors 
impede a CW gain output of the system. Therefore, the steady-state CARLs were proposed and have been 
experimentally demonstrated but, with a relative weak gain. In this paper, we report a result in order to give a 

more detailed theoretical analysis to a modified CARL system 10] which produces a CW laser output with a higher 
gain. In our previous scheme, the active cold atoms are trapped in a harmonic potential, which can be realized by 
loading the atoms into an optical lattice or an overall magnetic trap, instead of introducing a light molasses 
or adding a thermal collision/momentum damping mechanism to the atomic motion We suppose the trap to 
be strong in radial direction and loose in axial direction (cigar-shaped trap), which reducing our system effectively 
to one-dimension. Our modification by this configuration produces a very different laser gain profile illustrated by 
Fig. 2 10]. Fig. 2 displays a discrete spectral structure which is very similar to that of a nonlinear Jaynes-Cummings 



3x1 3 



CD 

1x10 3 









x=100 
























Ill 










. 1 . 






-28-24-20-16-12 -8 -4 4 8 12 16 20 24 28 



A 21 



FIG. 2: Long time gain profile G(t, A21) of the probe field in a trap with frequency of v — 2 and the time r = 100. The other 
parameters are the same as Fig.l for a convenient comparison. 

model of one atom in a high-Q cavity [T^ . When the pump-probe detuning satisfies a Raman condition A21 = nv, 
where n is an arbitrary integer and v is the effective axial trapping frequency, a distinct light amplification process 
occurs for certain n as displayed by Fig. 2. In this case, a self-organization of all the atomic oscillators along the 
trap axis is demonstrated by atomic distribution pattern and atomic orbit bunching in the phase space in this paper. 
The following contains our model with a more rigorous analysis of all the above interesting results in a framework of 
semiclassical method. 



II. SIMPLE MODEL AND GAIN SPECTRUM ANALYSIS 

The relevant system is supposed to be enclosed in a ring cavity 01 an d the operating medium is composed of 
a collection of N two-level atoms of mass m, moving along the axial direction z of a trap described by Vt rap — 
^m[u± 2 {x 2 + y 2 ) + v 2 z 2 \. The atoms are driven by a strong far off-detuned pump field propagating along z and 
probed by a counterpropagating weak probe field. If we set i>± 3> v z , the radial trapping mode will be hardly excited 
in the recoil process and a one-dimensional trap, V = ^mv z 2 z 2 , should be a promising approximation. With above 
simplification, the explicit form of the Hamiltonian of this system is 

N ( P 2 1 ^\ N ^ 

3=1 V m J 3=1 

(N N \ 

g{a\ dje- lk ^ + 9l a\ £ dje + ^ - H.c. , (1) 
3=1 3=1 / 

where the operators pj and Zj are the canonically conjugate momentum and center-of-mass position of the jth. atom, 
aj and a Z j denote the usual spin operators of the atomic internal degrees of freedom with a transition frequency of 
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ljq. The probe field and the pump field with the carrier frequencies toi and w 2 are denoted by the photon creation 



operators a \ and a 2 , respectively. The coupling constant between the atom and the ith optical field is given by <?j 
(i = 1,2). Following the conventional commutation relations of all above operators, the Heisenberg equations for the 
expectations of the relevant operators under the semiclassical approximation are 
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where the atomic polarization and population damping T as well as the probe field damping k have been phenomeno- 
logically added. The strong pump field A 2 is treated as a real undepleted field and removed from Eq.(J2J) as an 
external condition [2j. Some dimcnsionlcss variables or transformations have been introduced in above derivations, 
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^ ^e*"« (< = l,2),er i = <37 
A20 = {u 2 — uo)/(u r p) 1 v = v z /(uj r p) and t = uj r pt under the assumption k\ 

The scaling parameters u> r and p used above are two-photon recoil frequency u> r = 2hk 2 /m and the dimcnsionlcss 
CARL parameter p = I 9^. — J 1 respectively. It should be noted that Eq.@ are identical to those derived in the 

Ref.000 except for one single term —v 2 9 3 for the harmonic trapping force, and all the simulations throughout this 
paper will be based on these equations. 

In order to analyze the unique discrete gain profile in Fig. 2, we focus on the long-time evolution of A\ when the 
probe damping is weak so as to be enclosed in a high-finesse ring cavity [lflj . The gain quantity depicted in Fig. 2 

is defined by G(A2i,r) = ^|j^7ojpr^ with A\ (0) being the initial value of the weak probe field. By formally 

integrating Eq.(J2J), the probe field can be expressed as 
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, <j z j = -2 (a Z j), A 2 i = (ll> 2 - u>x)/(u> r p), 
k 2 = k and 171 « g 2 — g, for simplicity. 



At (r)=A 1 (0)e^ A 
and the probe gain 

G( K ,A 21 ,r) =e- 2KT 
where C(k — iA 2 i,r) is defined by 
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(5) 

with C(t) = i Y^j a j ( T ) e~ l9 j' r ' being defined as a collective atomic coherence parameter by Ref.@|. For an enough 
long time of evolution, provided r ^> 2n jv with v > 2p l^ lj42 l p^|. all the atoms in the trap will approximately perform 
a harmonic oscillation with different amplitudes 9jo(r) and phases 4>jo(r), described by 9j(r) « $jo( T ) cos [vt + (f>jo (r)]. 
Then e~ l9 ^ T ' in Eq.JSJ) will be expanded as 
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where the function J n [9jo(r)] is the n-th order Bessel function of the first kind. In above derivations, we have assumed 
that the amplitude 0jo(t') and phase (/>jo(i~') are slowly varying quantities after a long time evolution and can be 



4 



moved out from the integral of Eq.JBJ). If we extend the upper integral limit of Eq.© to infinity (r — > oo), the 
integral defines a Laplace transformation of a 3 - (r'), i.e., C [aj (r')] = / o °° dr'er., (r') • e l B - l ( A !i-™'")]'' ; on condition 
that \<jj (r') | doesn't increase much more than eF T . Generally <jj (r) is limited by Bloch sphere of two- level atom and 
can be approximately expanded by a Fourier series with v as a fundamental frequency 

oo 



°j ( r ) 



(7) 



where c„ is the Fourier amplitude of the n-th order harmonic components. Furthermore, if the pump-atom detuning 
satisfies A 2 o 3> T, the long-time average polarization (oj(r)) will be 
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where f2 = 2pA 2 is the Rabi frequency of the pump field and So can be determined by taking an adiabatic approx- 
imation of atomic internal freedom as done in |Tfij . If we only take the dc term of Oj (r) in Eq.JJJ) for a qualitative 
spectrum analysis, Eq.© reads 
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where Cq(t) = J2j=i Soe iej °( T ^ cos ^ j °( T ^ is the asymptotic coherence parameter which is determined by the secular 
average atomic polarization So and the instantaneous atomic distribution at r. The coefficient (1 — e~ KT ) /n of Ea. ll(J|) 
shows here that when the damping rate k is small enough and the evolution time r 3> 1/k, the intensity gain of the 
probe field will be extremely high in the neighborhood of A21 = nv, and thus results in a uniform discrete gain 
structure as shown in Fig. 2 simulated directly on Eq.©. Although the Raman gain condition will not be changed 
if the higher harmonic components of Eq.Q is included, except for a more complicated sum formula for Eq.Q, the 
height of the gain peak for a specific n in Fig. 2 can't be interpreted by above analytical method. As shown by Fig. 2, 
for example, the light amplification mostly for positive multiple of the trap frequency and approximately up to n = 10 
is still a problem. The further simulations only indicate that the light amplification for positive detuning is a feature 
owing to the red-detuned of the pump field from the atom frequency. However, this regular gain spectrum still can 
be developed as an efficient detection method to analyze the micro-trap potential formed by a local interference field. 



III. SELF-ORGANIZATION IN PHASE SPACE 



Under the Raman resonant condition A21 = nv, another ingredient that enhances the maximum gain of probe field 
is the synchronization of the atomic motion in the harmonic trap, which is indicated by the term of Cq(t)/Ai(0) in 
Eci. (|10fl and can be demonstrated through a bunching of the atomic trajectories in the phase space. To some extent, 
the present problem described by Eq. J2J) under the adiabatic approximation is similar to the famous Kuramoto model 



of ensembles of coupled oscillators Q , in which the complex order parameter i?(r)e i *( r ) = J2 jLi e l6j( - T ^ plays 
an equivalent role as that of the bunching parameter b in the traditional CARL system [la|7 For the complex 
order parameter, R(t) measures the phase coherence of all the oscillators and <I > (t) describes their average phase. A 
main result of Kuramoto model is: When the coupling intensity (in present model the coupling intensity is related 
to the Rabi frequency, pA 2 , of the pump field) exceeds a certain threshold, the incoherent oscillators are mutually 
synchronized and R(t) grows exponentially. In present case, when the resonant conditions A21 = n ■ v and A20 S> T 
are both satisfied, the long-time output gain of A± reduces to 
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where i?o(r)e l *°^ r ^ = J2j=i e t0jo ^ cos ^ jo ^ T ^ is the secular complex order parameter determined only by the atomic 
distribution in the phase space. Eg. 1)11(1 indicates that the field gain is explicitly dependent on the increase of collective 
atomic coherence, R(t). At an early stage of evolution, as the atoms incoherently scatter the pumping photon into 
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FIG. 3: Sample snapshot distributions of 1000 atoms in phase space (a) at r = (Initial distribution), (b) at r = 950, (c) at 
r = 1800 and (d) the corresponding time evolution of the probe intensity |Ai(t)| under the resonant condition with parameters 
A 2 i = 2u, v = 2, T = 1, k = 0.01, p = 3, A 2 = 2, and A 20 = -15. 

probe mode and their spatial distribution in the trap almost keeps uniform with R(t) ~ (Fig.3(a)), the probe field 
fluctuates around a small value. After a period of lethargy, a self-organization structure is established through the 
recoil motion modified by the trap and the atoms in phase space concentrate into several clusters (Fig. 3(b)), resulting 
in an exponential growth of the probe field. When this process is saturated (Fig.3(c)) by a counteraction of the 
spontaneous damping processes |l7|. the coherent parameter reaches its maximum value, or extremely R(t) — > 1 (all 
the atoms moving synchronously as one big atom), and a stable output of A\ will be set up. Fig. 3(d), which is directly 
simulated from Eq.Q, displays a typical amplification process of the probe intensity under this resonant condition 
and, finally, establishes a saturated output with an efficiency nearly to 18% (Here the efficiency is defined by the 
intensity rate between pumping field and the output probe field). According to a further Fourier analysis, the fast 
oscillation observed in Fig. 3(d) at saturation is induced by the harmonic motion of the atoms in the trap. 

During the lasing process, the self-organization structure related to the light amplification corresponds to a spon- 
taneous bunching of atomic trajectories in the phase space (Fig. 4), which, actually, is similar to the atomic density 
grating in the free space CARL system when b is large. In order to demonstrate the spatial features of this atomic 
self-organization behavior, the initial spatial distribution of the atoms are set uniform in a region of [0, A] and the 
initial atomic momentum follows a Gaussian distribution (showed by Fig. 3 (a)) with a thermal spread of width above 
recoil temperature (about 0.8phk) 01 . Before the increment of Ai, the atomic trajectories are almost ubiquitous in 
the phase space and, after a stable output of A\ is established, most atomic trajectories converge into the antinodc 
area (around 0.25A) and the other ones into the node area (around 0.5A) of a standing wave induced by an interference 
between the amplified probe field and the strong pump field, with its wavelength being defined by A = 47r/(fci + A^). 
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FIG. 4: The atomic trajectories in phase space during (a) 500 < r < 510 and (b) 1700 < r < 1710 with parameters v — 2, 
r = 1, k = 0.01, p = 3, A 2 = 2, A 20 = -15 and A 2 i = 4. 



The bunching positions of the most atomic trajectories into the antinode area are a clear result due to the red-detuned 
of the pump field A 2 o = — 15 < in Fig. 4. 

According to more simulations, we find that the amplification and the pattern formation of the atoms in phase 
space under the Raman resonant condition are robust to other parametric values provided that the gas medium is 
operated on a far-off red-detuned region and the trap is strong enough [Til ] . Following above conditions the change 
of collective coupling constant p in the resonant gain region will not qualitatively change the self-organization effect 
and the long time high CW gain output of this model. 

IV. DISCUSSIONS 

The experimental realization of above CW CARL system can be carried out by loading the laser-cooled atoms into an 
overall magnetic trap or into a far-off detuned optical lattice. Thereafter, the trap is displaced along its axis for about 
one light wavelength A to trigger the atomic motion, and then the atoms are pumped by a strong red-detuned field. 
A gain of a seed probe field in the backward direction under the resonant conditions A21 = 0, ±z/, ±2i/, ±3z/ • • • will be 
detected then. This process will be enhanced further by a collective behavior of the atoms due to the synchronization 
of their recoil motion under Raman resonant condition. The strong pump field feeds the energy to the system (excite 
the trapping oscillating mode by atomic recoil motion) and the dipole force, generated by the interference of the probe 
and pump light, bunches the atoms into several groups to adjust them synchronically extracting kinetic energy from 
the trap mode into the probe field Q . This scheme will also present a clue to introduce a new configuration of free 
electron lasers (FEL) where the high speed electrons are injected into a (or an array of) strong quantum trap (e.g. 
quantum dots). By stimulating the electrons' motion through electromagnetic field on the electrode with appropriate 
frequency, the high gain field with certain frequency will be generated when the electrons oscillate collectively in (or 
through) these quantum traps. 

In conclusion, we have investigated the light amplification in a CARL with the atomic center-of-mass motion 
influenced by an external harmonic potential. Under appropriate parametric condition that the pump field is red- 
detuned from the atomic transition frequency, our results show that CARL confined by a robust trap can give rise to 
a CW laser with a higher lasing gain than that of the CARL in the free space. The optical gain profile as a function 
of the pump-probe detuning displays a discrete structure and another kind of atomic self-organization related to the 
light gain is demonstrated in the phase space. This work suggests a new way to promote a CW CARL device with a 
higher gain and, perhaps, analogically implies a new configuration of FEL. 
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